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ABSTRACT 

It is well-known that galaxy environment has a fundamental effect in shaping its properties. We 
study the environmental effects on galaxy evolution, with an emphasis on the environment defined as 
the local number density of galaxies. The density field is estimated with different estimators (weighted 
adaptive kernel smoothing, 10 th and 5 th nearest neighbors, Voronoi and Delaunay tessellation) for a 
K s <24 sample of ~190,000 galaxies in the COSMOS field at 0.1<z<3.1. The performance of each 
estimator is evaluated with extensive simulations. We show that overall, there is a good agreement 
between the estimated density fields using different methods over ^2 dex in overdensity values. How¬ 
ever, our simulations show that adaptive kernel and Voronoi tessellation outperform other methods. 

Using the Voronoi tessellation method, we assign surface densities to a mass complete sample of qui¬ 
escent and star-forming galaxies out to z~3. We show that at a fixed stellar mass, the median color 
of quiescent galaxies does not depend on their host environment out to z~3. We find that the number 
and stellar mass density of massive (>10 n Mo) star-forming galaxies have not significantly changed 
since z~3, regardless of their environment. However, for massive quiescent systems at lower redshifts 
(z<1.3), we find a significant evolution in the number and stellar mass densities in denser environ¬ 
ments compared to lower density regions. Our results suggest that the relation between stellar mass 
and local density is more fundamental than the color-density relation and that environment plays a 
significant role in quenching star formation activity in galaxies at z<l. 


1. INTRODUCTION 

It is known that the environment in which galaxies 
reside plays a fundamental role in their evolution. In the 
local universe, denser environments are dominated by 
red, passive, early-type galaxies whereas less dense re¬ 
gions are preferen tially populated by blue, star-forming , 
late-type systems (iDresslerll 1 98(1 [Kauffmann et al.ll2d04 
Baloeh et all 12004 iBaldrv et al.l 120061: iBamford et al.l 


a alo gn et al.l !2UUlt IDali 

2009( IPeng et all 120101b 


These environmental trends 


still hold at higher redshifts (Capak et al.l 2007a; 

Cooper et al. 

2007 

Peng et al.l 20~ 

L0l Sobral et al. 2011: 

Muzzin et al. 

2012 

: Scoville et al. 

2013 ), although they 


usually tend to weaken with increasing redshift. While 
it is evident that almost any observable property of a 
galaxy demonstrates some association with environment, 
there is still a question that needs to be addressed first. 
What do we really mean by “environment” ? 

Recent advances in numerical simulations (Millennium , 
iSoringel et al.1 [20051 Illustris, iVogelsberger et ahl 120141 
combined with extensive spectrosco pic observations 
of local galaxies (SPSS, lYork et ahl I200C ; 2d FGRS, 
iColless et al.1 12004 GAMA. iDriver et al.l 201 If) have 


revealed that the u niverse has a web -like pattern, i.e.; 
the “cosmic web” (|Bond et alJ fl99 6) , containing dark 
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matter, gas and luminous galaxies. Galaxies in the 
cosmic web are organized in a network containing dense 
clusters, sparsely populated voids, planar walls, and 
thread-like filamentary structures linking overdense re¬ 
gions. Therefore, the most natural approach in defining 
the environment of a galaxy is to locate it within the 
cosmic web. However, the complexity and lack of a fully 
objective method in identifying the major components of 
the cosmic web often limit the environmental studies of 
galaxies within the comic web to numerical simulations 
or large spectroscopic surveys in the local universe. 

The conventional method of defining the environment 
as two extreme regions in the density distribution of 
galaxies, i.e. galaxy cluster and the general field, does 
not usually account for the full dynamical range of 
the density field. For example, there are intermedi¬ 
ate environments such as galaxy groups, outskirts of 


clusters and filaments which are ec 

ually important 

(Fadda et al. 20081: 

Porter et al. 2008t 

Tran et al. 20091; 

Geach et al.1 201 ll 

Coppin et al. 

20121; Mahaian et al. 

2012: Pintos-Castro et al. 2013t 

Darvish et al. 2014). 


Moreover, the selection of galaxies to whether they 
belong to the cluster environment or the general field is 
somewhat subjective. 

The detection of gal axy clusters through their dif¬ 
fuse X- ray emission 

19981: iBOhringer et al _ . . . . .. 

Bbhringer et al.l 12004 iMehrtens et al.l l2 012fb t he red 
sequence a nd i ts variants dGladders fe Yeel l2000l: 
Goto et all l2002t Gladders fe Yeel 120051 : iMiller et al 


axy clusters tnrougn tneir ait- 
dGioia et~ah 1199( 4 lEbeling et, al l 

1.1 1200(4 lEbeling et all 120011: 


200(L Koester et al.l 120071 : lHao et al l 1 201 Ol lk bizziii et al. 

2013a ), weak gravitation a l lensi ng ([Miyazaki et al.ll2007t 
Gavazzi fc Soucaiil 120071 : iDietrich et al.1 12 0071 a nd the 
Sunyaev-Zel’dovich effect dPlanck Collaboration et all 




























































































2 


201 It iWilliamson et al.1 120111 : iMarriage et al.1 120111 : 
Reichardt et al .1120131) adds more assumptions and limi¬ 
tations in the environmental studies. For example, weak 
lensing, hot X-ray emission and the Sunyaev-Zel’dovich 
effects are more sensitive to virialized, massive galaxy 
clusters. The red sequence technique, which relies 
on observations in only two filters, is very successful 
and economically efficient in detecting a large sample 
of galaxy clusters that have a tight red sequence of 
quiescent galaxies. However, this technique assumes the 
existence of a tight red sequence for clusters, requires 
modelling the red sequence and uses only the quiescent 
systems as a proxy of their host cluster environment. 
Moreover, the quiescent galaxies become less abundant 
at higher redsliifts which makes the cluster detection 
techniques based on the red sequence (or galaxy color) 
challenging at high redshift. All these methods may 
produce biased environments, which may lead to misin¬ 
terpretation of the evolutionary history of galaxies as a 
function of their host environment. 

Another approach in defining environment is to use 
the local number density (usually surface density) of 
galaxies as a proxy to their host region. Among the 
most common density estimators used in the literature 
are the nearest neighbor (usually 5 th or 10 th NN) and 
count-in-cell (CC) methods. However, although simple 
to implement and computationally fast, the performance 
of these methods strongly depends on the depth of the 
data, the number of neighbors considered in the analysis 
(in NN method) and the size of the cell (in CC method). 
A small value of N (a small cell size for CC method) 
results in a spiky density field which makes it vulnerable 
to unrealistic density values due to Poisson noise and 
random clustering of spatially uncorrelated galaxies. A 
large value of N (a large cell size for CC method) is 
prone to underestimation of the surface density and 
over-smoothing the details of galaxy distribution. Also, 
for the nearest neighbor method, the sum of the area 
(volume) assigned to each galaxy is not equal to the 
total area (volume) of any survey and it has also been 
shown that its integral over all area (volume) diverges. 
More importantly, it is still not clear whether the 
environmental effects depend on the physical scales 
at which the env ir onme nt is estimated. F or exa mple, 
Kauffmann et al.l (12004H . iBlanton et al.1 (1200611 and 
Cucciati et al.l (l201(lh find that the effect of environment 
on the star formation history of galaxies is only effect ive 
at the small scales (<1 Mpc) whereas [Balogh et al.l 
(j2004f) find that the equivalent width of Ha (a measure 
of the specific star formation rate in galaxies) is a 
function of environment measured on scales of ~1 and 
~5 Mpc, independent of each other. Scale-dependent 
density estimators such as the NN and CC algorithms 
are not able to directly address this issue. 

In this work, we perform a comprehensive analysis of the 
(surface) density field using the weighted versions of the 
5 th NN, 10 th NN, adaptive kernel smoothing, Voronoi 
tessellation and Delaunay triangulation methods. This 
is done using a K. 9 -b and selected sample of galaxies in 
the COSMOS field (jScoville et al.l 12007b ). The large 
size of the COSMOS (~2 deg 2 ), together with the 
high accuracy of the photometric redshifts in this field, 
enable us to delineate the Large Scale Structure with 


great accuracy out to z~3. In this analysis, we use the 
full photo-z probability distribution function (PDF) of 
individual galaxies to significantly reduce the projection 
effect. The performance of each method is checked with 
extensive realistic and Monte-Carlo simulations. We 
then apply the algorithm to a mass complete sample of 
galaxies to investigate the role of density-based defini¬ 
tion of environment on the rest-frame color evolution of 
quiescent galaxies, as well as the evolution of comoving 
number and mass density of massive systems out to 
z~3. 

The format of this paper is as follows. In Section 
O we briefly review the data used to estimate the 
density field. Section [3] outlines the algorithm used 
to determine the density field, followed by the surface 
density estimation methods in Section [4] Section [5] deals 
with the simulations used to check the performance 
of different density estimators. Comparisons between 
different methods are given in Section [6] In Section [71 
we study the color evolution of quiescent galaxies, as 
well as number and mass density of massive systems as 
a function of environment, using a stellar mass complete 
sample from the COSMOS field. We give a summary of 
this work in Section [8] 

Throughout this work, we assume a flat concordance 
ACDM cosmology with Ho=70 kms -1 Mpc -1 , S2 m =0.3 
and Da=0.7. All magnitudes are expressed in the AB 
system and stellar masses are given assuming a Cliabrier 
IMF. 

2. DATA & SAMPLE SELECTION 

In order to examine different density estimation tech¬ 
niqu es, we select an area of ~1.8 deg 2 in the COSMOS 
field dScoville et al.ll2007bl : lCanak et al. 11200711 ) . Here, we 
use the COSMOS UltraVISTA K s -band selected photo¬ 
metric redshift (p hoto-z) catalog dMcCracken et al.ll2012l : 
lllbert et "all 1201311 . The K s -band selection is equivalent 
to a stellar mass-selected sample and enables us to re¬ 
liably detect galaxies (especially quiescent systems) at 
intermediate to high redshifts. This catalog consists of 
ground and space based photometric data in 30 bands, 
spanning the range UV to rnid-IR. Using this extensive 
data set, accurate photometric redshifts are measured 
for galaxies out to z~3. Stars, AGNs and X-ray sources 
were identified and removed from the sample. The final 
selection criteria used to define the sample are: 

• An area of 149.3< 0:2000 <150.8 and 1.6< 

J 2000 <2.8 corresponding to ~1.8 deg 2 . The area 
is selected to be large enough for the effect of LSS 
to be discernible. 

• Redshift within the range 0.05<z<3.2 to allow pho¬ 
tometric redshifts with high degree of accuracy. 

• Galaxies brighter than K, <24, since the photo-z 
uncertainties increase significantly at K s >24. 

The total sample contains 191151 galaxies. 

3. DETERMINATION OF THE LARGE SCALE 
DISTRIBUTION OF GALAXIES 

In any study of the LSS of galaxies, we need the an¬ 
gular position (a, S) of galaxies and a measure of their 
radial distance (using spectroscopic or photometric red¬ 
shifts) to estimate the density field. Here, we use the full 
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Fig. 1. — The photo-z uncertainty of galaxies used in this study. 
The blue solid line shows the median of the photo-z uncertainty 
(Az med) at each redshift. Az for each galaxy is calculated as the 
(lower &; higher) 68% confidence interval of its photo-z PDF. Note 
the accuracy of photo-z values. The median of the photo-z uncer¬ 
tainty is Az me d <0.01 out to z~l, reaching Az mec j ~0.1 at z~2. 

photo-z probability distribution function (PDF) of each 
individual galaxy to significantly reduce the projection 
effect. In summary, the density field estimation involves 
the following steps: 

1. First, we construct a series of overlapping redshift 
slices (z-slices). The widths of these slices are ob¬ 
tained by the photo-z accuracy of the underlying 
sample of galaxies at each redshift. 

2. For each z-slice, we select a subset of galaxies be¬ 
longing to that slice. The galaxies in the subset 
are selected in such a way that the median of their 
photo-z PDF is located within that z-slice. For a 
given subset, we associate a weight to each galaxy. 
This weight shows the likelihood of each galaxy be¬ 
longing to the z-slice of our interest. The weight is 
determined by the photo-z PDF of each galaxy and 
the boundaries of the z-slice. 

3. At each z-slice and for its subset of galaxies, we 
then compute the weighted surface density field us¬ 
ing one of the estimation methods described in Sec¬ 
tion U 



requires a homogeneous set of multi-waveband photomet¬ 
ric data. Thanks to the wealth of the photometric data 
available in the COSMOS field, we have been able to 
obtain highly accurate photometr ic redshifts for a large 
population of galaxies in this field ( Mobasher et al.ll2007t 
lllbert et "all I2009LI2013H . Figure [T] shows the median of 
photo-z uncertainty (Az me£ z) as a function of (photomet¬ 
ric) redshift for the sample of galaxies used in this study. 
Az for each galaxy is calculated as the (lower & higher) 
68% confidence interval of its photo-z PDF. The median 
of the photo-z uncertainty is Az me d <0.01 out to z~l, 
reaching Az me d ~0.1 at z~2. Studying the 3D density es¬ 
timations solely based on the photo-z of galaxies, without 
taking their photo-z uncertainties into account, may re¬ 
sult in an erroneous evaluation of the density field. This 
is mainly due to the fact that the physical lengths that 
correspond to even the smallest photo-z uncertainties are 
larger than the typical sizes of most physical structures 
such as galaxy clusters and groups (e.g. at z~l, the 
photo-z uncertainty of Az^O.Ol is equivalent to the co- 
rnoving radial length of ^24 Mpc). Here, we limit our 
analysis to the 2D surface density estimations by con¬ 
structing a series of narrow slices in the redshift space. 
This is done by defining a series of overlapping redshift 
slices (z-slices). The width of each slice is set by the 
photometric redshift accuracy of the data at any given 
redshift. The slices overlap to allow proper contribution 
from galaxies which reside close to the boundaries of each 
slice. In practice, we make sure that approximately half 
of each z-slice is trespassed into its adjacent z-slices. The 
width of each z-slice (<5z) at redshift z is defined as twice 
the median of the photo-z uncertainty (Az me d) at that 
redshift: 

Sz — 2A Zmed (1) 


Figure |T] shows how Az me d changes with redshift. In this 
study, a total number of 133 overlapping z-slices are used 
over the redshift range 0.1<z<3.2. 

Choosing a much narrower z-slice results in smaller num¬ 
ber statistics and decreased completeness of the under¬ 
lying sample of galaxies. This subsequently leads to an 
underestimation of the surface density field. On the other 
side, broadening the z-slice increases the risk of contam¬ 
ination from foreground and background galaxies as well 
as a possible overestimation of the density field. Thus, 
our choice is a compromise between both limits. 


In the following, we give a detailed description of points 
1 & 2. Using this information, we then apply different 
density estimators to the data in Section [4] 

3.1. Redshift Slicing 

In order to construct the density field of galaxies, 
one needs a large, homogeneous and unbiased sample 
of galaxies with very accurate redshifts. However, there 
are serious complications in constructing such samples. 
Firstly, building a homogeneous sample is difficult. Sec¬ 
ondly, relying on a purely spectroscopic sample would 
bias the study as it would only concentrate on bright 
objects and mostly those with emission lines (i.e., star¬ 
forming galaxies which are often less clustered compared 
to quiescent systems). Spectroscopic samples also do 
not usually target enough number of galaxies in dense 
regions. Thirdly, finding accurate photometric redshifts 


3.2. Photo-z Probability Weights 

To estimate the surface density field, we use the full 
information of the photo-z PDF for individual galax¬ 
ies. First, for each z-slice, we select a subset of galaxies. 
Galaxies in each subset are selected in such a way that 
their median photo-z PDF is located within the bound¬ 
aries of that z-slice. We assign a weight (w*) to each 
galaxy in this subset based upon its photo-z PDF. This 
probability weight is calculated by measuring what per¬ 
centage of the photo-z PDF of each galaxy is contained 
within the boundaries of the considered z-slice. This 
weight represents the likelihood that the galaxy in ques¬ 
tion resides within that z-slice. In practice, it is done by 
intersecting the z-slice of our interest with the photo-z 
PDF of each galaxy in the subset. In order to enhance 
the computational performance, we assume a Gaussian 
photo-z PDF for galaxies except for those that have a 
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Table 1: Overdensity values evaluated using the adaptive kernel 
estimator for sample galaxies selected in Section [2] 


ID 

<*2000 

deg 

<52000 

deg 

photo-z 

K s 

mag 

log(£/Emedian ) kernel 

1 

150.02025 

1.67951 

0.0532 

22.69 

-0.5811 

2 

150.39144 

1.68383 

0.0506 

23.30 

-0.6667 

3 

149.58068 

1.68829 

0.0512 

23.71 

-0.4306 


TABLE 1 

Notes: Table 1 is published in its entirety in the 
electronic edition OF the Astrophysical Journal. A PORTION IS 
SHOWN HERE FOR GUIDANCE REGARDING ITS FORM AND CONTENT. 


second photo-z peak with P>5% in their PDF. For the 
latter, the full photo-z PDF is used in estimating the 
weights. Eventually, at each redshift and for each z-slice, 
we have a subset of galaxies with associated weights that 
are later used as the input for the weighted surface den¬ 
sity estimation methods. Figure [2] (a) shows the distri¬ 
bution of galaxies for a z-slice centered at z=0.92 in the 
COSMOS. For demonstration, the size of each point is 
selected to be proportional to the galaxy weight. This 
subset of galaxies with their associated weights is used 
to estimate the weighted surface density field at z=0.92. 

We note that the weight introduced in this section 
can be generalized to other galaxy properties such as 
stellar mass, luminosity and color, depending on the as¬ 
trophysics of our problem. 

4. SURFACE DENSITY ESTIMATION: THE METHODS 

Using the information obtained in Section [3l we now 
construct the density field for galaxies in the COSMOS, 
using weighted versions of four independent methods: 
adaptive kernel, nearest neighbor, Voronoi tessellation 
and Delaunay triangulation (Sections 14.II to 14.411 . We 
then evaluate their performance with simulations in Sec¬ 
tion [5] and compare them with each other in Section [G] 


4.1. Weighted Adaptive Kernel Estimator 

In this method, we use an iterative procedure to com¬ 
pute the surface density field dDarvish et al.ll2014l ). First, 
we estimate the surface density associated with the *th 
galaxy in a given z-slice, E^, by summing over all the 
weighted fixed kernels placed on the positions of galax¬ 
ies, j, where i ^ j: 

1 N 

±i = —^ -J] wj K (r i; rj , h) (2) 

1 W 3 j — I 

jA i 


where N is the number of galaxies in the subset, 
iF(r;,rj,h) is the fixed kernel, r; is the position of the 
galaxy for which the initial estimate of surface density is 
measured and rj is the position of the rest of the galaxies. 
The width of the kernel is expressed by the parameter, 
h, which is a proxy for the degree of smoothing. For the 
first estimate of the density, this is taken to be fixed. For 
the kernel smoothing function, we use a 2D symmetric 
Gaussian defined as: 


/C( r i,rj,h) 


1 

27rh 2 



( 3 ) 


A large kernel width (h) results in over-smoothing of the 
density field which tends to wash out real features while 


a small value tends to break up regions into smaller un¬ 
correlated substructures. Here, we use a fixed physical 
length of h=0.5 Mpc which corresponds to the typical 
value of R? nn for X-ray clusters fc groups in the COS - 
MOS field IlFinoguenov et alll2007t iGeorge et all 1201 111 . 
However, a constant value of h for the whole field has 
the problem that it underestimates the surface density 
in crowded regions while overestimates in sparsely pop¬ 
ulated areas. To overcome this problem, we introduce 
adaptive smoothing width, hq, which is a measure of 
the local surface density associated with each galaxy, E^. 
This is defined as hj=hxAi, where A i is a parameter that 
is inversely proportional to the square root of the surface 
density asso ciated w ith the i th. galaxy, at the position of 
that galaxy (jSilvermanlll986 ): 

Ai = (G/ECd)) 0 ' 5 (4) 

Where G is the geometric mean of all the E(ri) values. 
Having the adaptive kernel, we now calculate the surface 
density field, E(r), on each location on a fine 2D grid, 
r=(x,y) as: 

i N 

S(r) = — N - Y! w i K (u r i> h i) (5) 

L i= i wi i=1 

The surface density field is evaluated on a fine grid with 
a grid size (resolution) of 50 Kpc at each redshift. We 
note that the surface density field estimated through this 
method is almost independent of the type of the kernel 
function. We examined several other standard kernels in¬ 
cluding exponentially decaying, Epanechnikov (K(r)oc(l- 
r 2 )) and cosine arch (K(r)occos(^r)) and did not find any 
significant difference in the final results. Throughout this 
work, we define overdensity as: 

1 + 5 = ^— ( 6 ) 

where E median is the median of the surface density field 
at each redshift. Figure [2] (b) shows an example of the 
surface density field for a z-slice in the COSMOS field 
centered at z=0.92, using the weighted adaptive kernel 
density estimator. Note the wide range of overdensi¬ 
ties and the variety of environments, including dense 
regions linked together through thread-like filamentary 
structures. 

In table 1, we have provided an example of the overden¬ 
sity values estimated with the adaptive kernel smoothing 
method for our sample galaxies (Section [2]). K s mag¬ 
nitudes and photo-z estim ates are extracted from pub- 
licly available cat alogs of (IMcCracken et all 120121) and 
(|Ilbert et all l2013h . respectively. We recommend using 
galaxies that are far from the edge of the field and/or 
the masked regions. A catalog of overdensity values with 
other density estimation methods is also available upon 
request. 

4.2. Weighted k-Nearest Neighbors(k-NN) estimator 

In a regular (non-weighted) k-NN method, the inverse 
of the area containing the kth nearest neighbors to each 
galaxy is used as a proxy for local surface density at 
the position of that galaxy. Since galaxies in our sample 
are weighted, we incorporate the role of weights in the 
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Fig. 2.— Comparison of surface density fields in the COSMOS at z=0.92, estimated with different methods, (a) Distribution of galaxies 
for a z-slice located at z=0.92. The size of each point is proportional to the galaxy weight, (b) Surface density field estimated using 
the weighted adaptive kernel smoothing method for the distribution of galaxies at z=0.92. (c) Surface density field estimated using the 
weighted 10 th NN method (d) 5 th NN (e) Voronoi tessellation and (f) Delaunay triangulation methods. 


density estimation. The surface density at the position 
of the ith galaxy, E(ri), is estimated as: 


E(ri) 


EjLijwij 

Ej k =i Wijdfj 


(7) 


Where r; is the position of itli galaxy, w ij is the photo-z 
probability weight (Section 13.21) associated with the jth 
nearest neighbor to the galaxy located at r, and d,j is 


the distance between them (the distance between the 
galaxy positioned at r, and its jth neighbor). Using 
the distance to the fifth (k=5) or tenth (k=10) nearest 
neighbor is very common in litera t ure (s ee for example, 
iCooper et al.l (120051 ) : ISobral et all (| 2011 f )b In this work, 
we use the weighted versions of both the fifth (NN 5 ) and 
tenth (NN 10 ) nearest neighbor estimators. Figures [2] (c) 
and (d) show examples of the surface density field at 
z=0.92, estimated using NNi 0 and NN 5 , respectively. 
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The NN 5 method gives a larger density value in very 
dense regions compared to NN 10 , as expected. The over¬ 
all traces of NN5 are seen in NN5 plot. However, the 
NNs-based density field looks spikier and clumpier due to 
the smaller physical lengths that the NN 5 method spans. 

4.3. Weighted Voronoi Tessellation Estimator 

In a simple Voronoi tessellation method, we divide 
the 2D z-slice plane into a number of regions (Voronoi 
cells) assigned to each galaxy located at 13 . The 
Voronoi cell of each galaxy is defined as all points in 
the z-slice plane that are closer to that galaxy than 
to any other galaxy. It is acquired from the intersec¬ 
tion of half-planes. Based on this definition, in more 
crowded, denser regions, the Voronoi cells of galax¬ 
ies incline towards smaller values, therefore, we can 
use the inverse of the area of the Voronoi cell of 
each galaxy as a measu re of t he local dens ity at the 
posit ion of that galax y (lEbelin g & Wiedenmannl 11993! 

[Bernardeau fc van de Wevgaertlll996[ ) This is given by: 

£(r0 = ^ (8) 

Where E(ri) is the surface density at the position of the 
ith galaxy and A* is the area of its Voronoi cell. 

In order to assimilate the role of galaxy weights, we tried 
two different approaches. In the first approach, we mod¬ 
ified the metric used for the definition of distance be¬ 
tween points in the z-slice plane and galaxies. This was 
done by dividing the Euclidean metric by the weight of 
each galaxy. Implementing this approach to the data was 
not successful and resulted in unrealistically high density 
values for the majority of galaxies (with small weights) 
located in sparsely populated regions. Thus, we alterna¬ 
tively use a Monte-Carlo acceptance-rejection process in 
order to take the role of weights into consideration. The 
steps are as follows: 

1. For each galaxy with its corresponding weight w,, 
we generate a random number R, between the min¬ 
imum and maximum weight values in each z-slice. 

2. If Wj > R,:, we accept the galaxy (with its asso¬ 
ciated weight Wj) in the density estimation proce¬ 
dure. 

3. We estimate the surface density for the ac¬ 
cepted galaxies using a simple Voronoi tessellation 
method. This surface density is evaluated at the 
position of accepted galaxies only. 

4. We interpolate^ the estimated densities into the 
points of a grid (r) (the grid resolution is 50 Kpc 
at each redshift) constructed on the z-slice plane. 
This gives us a Monte-Carlo estimated surface den¬ 
sity field E(r). 

5. We repeat the above procedure N times and take 
the mean of all the Monte-Carlo density fields as 
the actual density field E(r): 

S(r) = 

’ m=l 

1 We use natural neighbor interpolation developed by Sibson. 


Finally, the local surface density of each galaxy is then es¬ 
timated as that of its closest point in the grid constructed 
over the z-slice plane. To save computational time, we 
use N=10 in step 5. Figure [2] (e) shows an example of 
the surface density field at z=0.92, estimated with the 
weighted Voronoi tessellation algorithm. Note the large 
dynamical range of overdensities. Unlike the nearest 
neighbor, Voronoi tessellation is scale-independent and 
is able to span a wide range of physical lengths. Also, it 
does not make any assumptions about the geometry and 
morphology of the structures in the density field. This 
characteristic makes it superior to adaptive kernel and 
nearest neighbor methods. 

4.4. Weighted Delaunay Triangulation Estimator 

This method relies on segmenting the z-slice plane into 
triangles whose vertices are defined by the position of 
galaxies in the z-slice plane. For each triangle, these 
three vertices (position of galaxies) are selected such that 
their circumcircle does not encompass any other galaxy. 
In this method, each galaxy is eventually surrounded by a 
series of neighboring triangles w hose overa ll area tend to 
be smaller in denser regi ons (iSchaap fe van de Wevgaertl 
120001; ' Platen et al.ll201lh . For the ith galaxy surrounded 
by to neighboring triangles, the estimated surface density 
is expressed as the inverse of the sum of the areas of its 
neighboring triangles: 

E( ri ) = 1 (10) 

2_m=l a n 

Where E(r;) is the surface density at the position of the 
ith galaxy and & n is the area of the nth triangle neigh¬ 
boring the ith galaxy. 

In order to take the weight of galaxies into account, we 
utilize a Monte-Carlo acceptance-rejection approach ex¬ 
plained in Section l4~3l Figure [2] (f) shows an example of 
the surface density field at z=0.92, estimated with the 
weighted Delaunay method. When compared with other 
estimation methods, this method overestimates the den¬ 
sities in very dense regions and the resulting density field 
is clumpier and contains much more substructures. 

5. SIMULATIONS 

In order to evaluate the performance of each surface 
density estimation method, we run two sets of simula¬ 
tions. The details are given in the following subsections. 

5.1. Simulation 1 

In the first set of simulations, we apply the surface 
density estimators to a sample of galaxies randomly 
draw n from some previous ly known surface density pro¬ 
files (|Scoville et a.l.l [201)7 all . Since the simulated surface 
density values are known, we can directly compare them 
with the estimated surface densities that are predicted 
by each of our density estimators. Here, we make 30 dif¬ 
ferent simulated structures. These structures are placed 
on an area similar to that of the COSMOS field. The 
simulated structures and the field contain 5000 galaxies 
at z = 1 , similar to the total number of observed galaxies 
at that redshift in the COSMOS field. For simplicity, 
we assume all to have the same weight (w=l). 2222 of 
these galaxies are randomly distributed on an area cover¬ 
ing 1.5x1.2 deg 2 . The rest of them are drawn randomly 
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Fig. 3.— (a) Overdensity map for simulated structures. Here, we make 30 different simulated structures placed on a constant background. 
They are assumed to be located at z=l and cover an area similar to that of the COSMOS field. The structures have Gaussian profiles with a 
variety of sizes (0.1-2 Mpc) and the number of galaxies in these structures is in the range 3-300. (b) Distribution of galaxies randomly drawn 
from the structures. The simulated structures and the field contain 5000 galaxies. Almost half of them (2778) belong to the structures 
while the rest are randomly distributed on an area covering 1.5x1.2 deg 2 , which serve like the field, (c) Predicted overdensity maps using 
the adaptive kernel, (d) 10 th NN, (e) 5 th NN, (f) Voronoi tessellation and (g) Delaunay triangulation. All these methods perform relatively 
well when compared to the expected density field, with adaptive kernel and Voronoi tessellation performing relatively better than other 
estimators. 




































from a series of Gaussian profiles. These Gaussian struc¬ 
tures have a variety of sizes (0.1-2 Mpc) and the number 
of galaxies in these structures is in the range 3-300. The 
properties of these structures are shown in table 2. We 
apply the density field estimation methods explained in 
Section 0] to the simulated galaxies to explore how well 
they can predict the input density field. Figure 0] shows 
the density field of the simulated structures, the distri¬ 
bution of galaxies randomly drawn from the structures 
and the predicted density fields using different estima¬ 
tion methods. All these methods perform relatively well 
when compared to the expected density field. We stress 
that some of the difference between the expected and pre¬ 
dicted density fields is due to the shot noise, which be¬ 
comes specifically important when the number of galax¬ 
ies in simulated structures is small and the width of the 
profile is large. In order to quantify the performance of 
different density estimation methods, we define the mean 
squared error (MSE) as: 

1 N 

MSE=-Y J {K m -^i) 2 ( 11 ) 

V i=1 

Where N is the total number of generated galaxies, 
Af* m =log(l +8f lrn ) is the logarithm of the overdensity 
of the ith galaxy in the simulation and A* is a similar 
quantity, predicted by one of the density estimator meth¬ 
ods. A smaller value of MSE shows a higher similarity 
between the simulation and predicted density fields and 
is a measure of the performance of the estimation meth¬ 
ods in predicting the true value of the density field. The 
MSE values are tabulated in table 3. It is clear that the 
adaptive kernel, Voronoi tessellation and NN i0 methods 
outperform the NN 5 and Delaunay triangulation algo¬ 
rithms. 

We note that due to the catastrophic failure in esti¬ 
mating the photometric redshifts, the actual discrep¬ 
ancy between the simulated and predicted density fields 
can be higher.However, the catastrophic failure fraction 
(defined as the fraction of galaxies which satisfy |z p hot- 
z spec|/(l+z S pec)>0.15) for K s <24 galaxies in the COS¬ 
MOS field is relatively small and estimated to be <1% 
for b right so urces , reaching < 10 % for fainter galaxies 
filbert et al.ll2013lh We redo the simulation taking into 
account the catastrophic failure and estimate the new 
MSE statistics. We assume that 10% of simulated galax¬ 
ies are affected by the catastrophic failure (i.e., they do 
not actually belong to the simulated z-slice) and they 
are randomly distributed over the area containing the 
simulated structures (random distribution is a fair as¬ 
sumption since catastrophic failure mostly affects fainter 
galaxies which tend to be less clustered). This means 
that the actual background for the overdensity map of 
simulated structures should have a smaller value. With 
the new simulation, all density field estimators perform 
almost equally worse with catastrophic failure (as pre¬ 
dicted). The MSE statistics is increased by ~0.01-0.02 
(~3-6%) for all density estimators. In this work, we do 
not correct for the catastrophic failure but we highlight 
that it has a small effect and the trends are not affected 
much by it. 

5.2. Simulation 2 


Table 2: Properties of simulated structures. The structures have a 
Gaussian profile with a variety of sizes and number of galaxies. 


Structure 

£*2000 

deg 

^2000 

deg 

<7 

Mpc 

N 

1 

149.5 

1.8 

0.1 

300 

2 

149.7 

1.8 

0.2 

300 

3 

149.9 

1.8 

0.5 

300 

4 

150.1 

1.8 

1.0 

300 

5 

150.3 

1.8 

1.5 

300 

6 

150.5 

1.8 

2.0 

300 

7 

149.5 

2.0 

0.1 

100 

8 

149.7 

2.0 

0.2 

100 

9 

149.9 

2.0 

0.5 

100 

10 

150.1 

2.0 

1.0 

100 

11 

150.3 

2.0 

1.5 

100 

12 

150.5 

2.0 

2.0 

100 

13 

149.5 

2.2 

0.1 

50 

14 

149.7 

2.2 

0.2 

50 

15 

149.9 

2.2 

0.5 

50 

16 

150.1 

2.2 

1.0 

50 

17 

150.3 

2.2 

1.5 

50 

18 

150.5 

2.2 

2.0 

50 

19 

149.5 

2.4 

0.1 

10 

20 

149.7 

2.4 

0.2 

10 

21 

149.9 

2.4 

0.5 

10 

22 

150.1 

2.4 

1.0 

10 

23 

150.3 

2.4 

1.5 

10 

24 

150.5 

2.4 

2.0 

10 

25 

149.5 

2.6 

0.1 

3 

26 

149.7 

2.6 

0.2 

3 

27 

149.9 

2.6 

0.5 

3 

28 

150.1 

2.6 

1.0 

3 

29 

150.3 

2.6 

1.5 

3 

30 

150.5 

2.6 

2.0 

3 


In the previous simulations, although the “true” den¬ 
sity values are known (which make the comparison with 
the predicted densities feasible), the complex geometry 
of the real astronomical density fields and their diverse 
dynamical range are not fully considered; i.e, the simu¬ 
lated density profiles have simple mathematical shapes 
whereas the real astronomical data consist of an irregu¬ 
lar web of filaments, voids, walls and clusters with a va¬ 
riety of physical scales and geometries. Furthermore, the 
previous simulation greatly suffers from the shot noise, 
especially when the number of randomly drawn galaxies 
is small and/or the known density profiles have broad, 
close-to-the-field shapes. In order to consider the com¬ 
plexity of the cosmic web in the analysis, one can use 
the numerical simulations which produce mock galaxy 
distributions that resemble the real density maps. Since 
the “true” density field is unknown in these quasi-real 
simulations, we should rely on an estimator to have an 
initial guess of the density field. In the second set of 
simulations, we make a sample of galaxies that resemble 
the true distribution of galaxies in the COSMOS. This is 
performed using a Monte-Carlo acceptance-rejection ap¬ 
proach, a method similar to the one explained in Section 
14.41 taking the following steps: 

1. For each z-slice, we estimate the surface density 
field for the real data using one of the density es¬ 
timators explained in Section [4] For simplicity, we 
assume that all galaxies in the z-slice have the same 
weight (w=l). The estimated density field gives us 
a first order approximation of the shape of the true 
density field. 

2. For each z-slice, we randomly select a position, 
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Table 3: Comparison between the simulated density field (simulation 
1) and the one predicted by different estimators. We use the MSE 
measure for comparison. A smaller value of MSE indicates a better 


performance. 


Kernel 

NNio 

nn 5 

Voronoi 

Delaunay 

MSE 0.26 

0.28 

0.30 

0.27 

0.34 


Table 4: Comparison between the simulated density field (simulation 
2) and the one predicted by different estimators. We use the MSE 
measure for comparison. A smaller value of MSE indicates a better 
performance. 


Kernel 

NN 10 

NN 5 

Voronoi 

Delaunay 

MSE 0.003 

0.035 

0.062 

0.025 

0.158 


(x,y )random, in the area covering the data and as¬ 
sign a random number, E ra n dom, between the min¬ 
imum and the maximum density values in that z- 
slice to the point (x,y )random- We also report the 
density value at this random position from step 1 
and call it E. 

3. If E>E rQrldom , we accept the point (x,y) random as 
one of the points in the simulated dataset and con¬ 
sider E as its true density. 

4. For each z-slice, we repeat steps 1-3 until we make a 
simulated dataset with the same number of galaxies 
in the actual data set. 

Given the above steps, for each z-slice, we now have a 
sample of simulated galaxies with known surface densi¬ 
ties. The distribution of these galaxies resembles the ac¬ 
tual data in the COSMOS. We apply our surface density 
estimation methods to the simulated data with known 
densities and compare their expected and predicted den¬ 
sities. Although now, we consider the complexity of real 
galaxy distribution into account, these simulations rely 
on an initial estimator to determine the rough shape of 
the density field (i.e., step 1). We therefore use the adap¬ 
tive kernel density estimator in order to perform the first 
step of the simulation 0 . Table 4 shows the performance 
of different surface density estimators when compared 
with the new set of simulations. Here, we use the MSE 
measure introduced in Section O to compare the re¬ 
sults. According to the new simulations and the results 
of the MSE statistics, the adaptive kernel, Voronoi tes¬ 
sellation and NNio methods perform better in estimating 
the simulated density field. In both sets of simulations, 
the NN .5 and (especially) the Delaunay tessellation algo¬ 
rithms seem to fall behind compared to other estimators. 

6. COMPARISON 

Using two sets of simulations, we showed that all den¬ 
sity estimation methods perforin relatively well in esti¬ 
mating the density field. However, the adaptive kernel, 
Voronoi tessellation and NNi 0 are found to outperform 
the others. In this section, we compare different estima¬ 
tion methods together regardless of their performance. 
Due to the edge effect, the surface density of galaxies 

2 This results in a bias in favor of the performance of the adap¬ 
tive kernel density estimator when we compare it with simulations. 
However, we used other estimators as well to initiate the simula¬ 
tion and found that the adaptive kernel estimator has an overall 
relatively good performance in all these simulations. 



log(l + S)nnio 



-1 0 1 
log(l + S ) NN5 





Fig. 4. — Comparison between different estimation methods using 
176893 galaxies at 0.1<z<3.2. In order to minimize the edge effect, 
we limit the comparison to galaxies that are 1 Mpc away from 
the edge of the field and large masked regions. There is a good 
agreement between different estimation methods over ~2 orders of 
magnitude in overdensity values. However, when we compare the 
Delaunay triangulation method with the rest, it overestimates the 
density values in dense regions while underestimates the density in 
sparsely populated areas. 


located near the edge of the survey and masked areas 
are unrealistically underestimated. In order to take this 
into account, we limit the comparison to galaxies that 
are more than 1 Mpc (physical) away from the edge of 
survey and large masked areas. Figure U shows a compar¬ 
ison between different estimation methods using 176893 
galaxies at 0.1<z<3.2. There is a relatively good agree¬ 
ment between different estimation methods over ~2 or¬ 
ders of magnitude in overdensity values. However, when 
we compare the Delaunay triangulation method with the 
rest, it overestimates the density values in dense regions 
while underestimates the density in sparsely populated 
regions (field) when compared with adaptive kernel and 
NNio methods. 


















10 


z—0.1-3.2 



«2000 {deg) 


Fig. 5.— Projected overdensity map in the whole COSMOS field at 0.1<z<3.2, using the weighted Voronoi tessellation method. In 
constructing this, the overdensity maps from all the z-slices between z=0.1-3.2 are stacke d and normalized to the peak. We compare 
it with the stacked overd ensity map from Voronoi-based algorithm of IScoville et alj d2013ll . Contours are used to demonstrate the map 
f rom Scoville et alJ 12013 1. Contours are at levels 0.2 to 1 with 0.05 spacing between levels. There is a very good agreement between 
lScoville^mn ~ if20Tj ) and our work. We also find a relatively good a greement between the denser regions in our work and the position 
of X-ray clusters/groups (Finogucnov ct al. 2007; George ct al. 2011), as well as the projected mass map form weak lensing analysis of 
IMassev et all l200% in the COSMOS. 


Each surface density estimator method has its own ad¬ 
vantages and disadvantages as explored below: 

• Adaptive Kernel: Since we smooth the density field 
with a suitable kernel function, this method is less 
affected by the shot noise and possible random 
clustering of foreground and background sources. 
Calculation of the global kernel width (h) is mo¬ 
tivated based on some physical scales (e.g. the 
typical size of galaxy clusters and groups) which 
makes it suitable for practical observational situa¬ 
tions. We can easily add weights to the estimator 
by multiplying the kernel function with a proper 
weight. This method conserves the total number 
of galaxies; i.e, the integral of the surface density 
field over the whole area yields the total number of 
galaxies. However, the selection of the appropriate 
kernel size is a serious problem, as a small kernel 
size tends to overestimate the sparsely populated 
regions in the density field whereas a large kernel 
size washes out real features. This is partly over¬ 
come by adaptively changing the kernel width, at 
the expense of increasing the computational time 
by adding an extra step to initially estimate the 


density field at the position of galaxies. Since we 
estimate the surface density at the positions of a 
grid, an intermediate interpolation is required to 
assign the density values from the grid points to 
the position of galaxies. It is also computationally 
expensive in its adaptive format as it requires extra 
steps to determine the local adaptive kernel width. 

• k-NN: It is by far the easiest method to imple¬ 
ment and fastest to perform computationally. It 
estimates the surface density directly at the posi¬ 
tion of galaxies which does not require any inter¬ 
polation. Adding weights to this method is read¬ 
ily done. However, its performance strongly de¬ 
pends on the number of neighbors considered in the 
analysis (k). A small value of k results in a spiky 
density field which makes it vulnerable to unreal¬ 
istic density values due to Poisson noise and ran¬ 
dom clustering of spatially uncorrelated galaxies. 
A large value of k is prone to underestimation of 
the surface density and oversmoothing the details 
of galaxy distribution. The sum of the area as¬ 
signed to each galaxy is not equal to the total area 
of survey. It has also been shown that its integral 
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over all area diverges. 

• Voronoi Tessellation: This method covers a wide 
dynamical range in densities and is able to esti¬ 
mate densities in a broad range, from the dense core 
of clusters to sparse regions devoid of galaxies. It 
is a non-parametric and scale-independent method 
which makes no prior assumption about the shape 
of the density field. We directly evaluate the sur¬ 
face density at the position of galaxies. However, 
we cannot assign closed Voronoi areas to galaxies 
near the edge of the field. Adding weights to this 
method is not straightforward. Here, we used a 
Monte-Carlo method to take the role of weights 
into consideration. However, this comes at the ex¬ 
pense of a computationally expensive process by 
making several Monte-Carlo samples. Apart from 
its computational time, it is a robust estimator. 

• Delaunay Tessellation: In terms of advantages and 
disadvantages, it is very similar to the Voronoi tes¬ 
sellation method. However, we can assign delau- 
nay area to all galaxies in the sample even those 
that are located at the edge of the field. Also in 
this method, the total area assigned to galaxies sur¬ 
passes the area of the survey. Despite its similarity 
to the Voronoi tessellation, this method does not 
perform well compared to other estimation meth¬ 
ods (at least for the COSMOS field). 

Due to the good performance of the Voronoi tessellation 
method, its large dynamical range, scale independence 
and no prior assumption it makes about the morphology 
of the structures, we use it for the scientific analysis in the 
next sections. Figure [5] shows the projected overdensity 
map in the whole COSMOS field at 0.1<z<3.2, using the 
weighted Voronoi tessellation method. In constructing 
this, we stack the overdensity maps from all the z-slices 
between z=0.1-3.2 and normalize it to the peak of stacked 
overdensities. We compare it with the sta cked overden- 
sity m ap from Voronoi-based algorithm of lScoville et al.1 
(120131). Contours ar e used to demonstrate the map from 
IScoville et al.1 (120131) . Contours are at levels 0.2 to 1 with 
0.05 spacing between levels. I n terms of sample sele ction, 
our work is similar to that of lScovil le et al.l (120131). We 
see a very good agreement between IScoville et al.l (l2013h 
and our work. We also find a relatively good agreement 
between the denser regions in our work and the posi¬ 
tion of X-rav clus ters/groups (|Finoguenov et all 120071 : 
IGeoree et al.l 120 111) . When compared with the projected 
mass map form weak lensing analysis of iMassev et al.1 
(|2007fl in the COSMOS fie l d, our projected map agrees 
with that of IMassev et al.l (1200711 . We are able to re¬ 
cover all the massive structures in the weak lensing map. 
Our estimated density field is also consiste nt quantita¬ 
tively with the density field of IKovac et al.l (120101) who 
used zCOSMOS spectr oscopic data out to z~l, t he op¬ 
tical galaxy groups of dKnobel et al.1 12009112012H using 
zCOSMOS d ata set and the protocluster candidates of 
IChiang et al.l (120141) at 1.6<z<3.1 in the COSMOS field. 
In the following sections, we use the weighted Voronoi 
tessellation method to study the dependence of the ob¬ 
servable parameters on environment. 

7. DEPENDENCE OF THE PHYSICAL PROPERTIES OF 
GALAXIES ON THE ENVIRONMENT 
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Fig. 6.— Mass completeness limit for all the galaxies, along with 
that of the quiescent systems as a function of redshift. Mass com¬ 
pleteness limit is defined in such a way that only less than 5% of 
galaxies could be missed in the lower mass distribution of galaxies. 
Using the mass completeness limits for the quiescent systems, we 
define six mass complete samples at different redshifts, as shown 
with labels 1-6 here. 

Using the surface density field constructed in the COS¬ 
MOS field, we now study the dependence of the observ¬ 
able parameters of galaxies on their local density. We 
first need to define a mass complete sample of galaxies 
at different redshifts. In the following sections, we select 
the sample first and then use it to investigate the role of 
environment on properties of galaxies. 

7.1. Sample Selection 

In the following subsections, we first define a mass com¬ 
plete sample at different redshifts. We also explain the 
selection of quiescent/star-forming galaxies. 

7.1.1. Stellar Mass Complete Samples 

The magnitude cut K s <24 introduced in Section [2] 
results in a magnitude-limited sample whose mass com¬ 
pleteness is a function of redshift and stellar M/L ratio. 
In order to estimate the mass comp leteness o f o ur sam ple, 
we use the method explained in lPozzetti et al.1 (12010( 1 and 
lllbert etldl (120131) . First, for each galaxy, we calculate 
the limiting stellar mass (M/i m ); i.e., the mass it would 
have at its redshift, if its apparent magnitude were equal 
to the magnitude limit of the sample (K s =24). This is 
given by log(M ;im /M 0 )=log(M/M 0 )+O.4(K s -24), where 
M is the estimated stellar mass of the galaxy with ap¬ 
parent magnitude K s . At each redshift, the stellar mass 
completeness limit (M comp ) corresponds to the mass with 
95% of the galaxies having their Mu m below the stel¬ 
lar mass completeness limit. This guarantees that only 
less than 5% of galaxies could be missed in the lower 
mass regime of the stellar mass distribution of galaxies. 
The mass completeness limit also depends on the stellar 
M/L ratio and is higher for quiescent galaxies. Therefore, 
in constructing the mass complete samples at each red¬ 
shift, we make sure that we rely on the quiescent galax¬ 
ies to estimate the stellar completeness limit. This mini¬ 
mizes the loss of low-mass quiescent galaxies especially at 
higher redshifts. The selection of quiescent/star-forming 
galaxies is explained in Section 17.1.21 Figure |H] shows 
the stellar mass of galaxies as a function of redshift for 
magnitude-limited sample (K s <24) defined in Section[2j 
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Table 5: Properties of mass complete samples used in this study 


Sample Number 

Redshift Range 

Mass Completeness limit 
log(Aig) 

Number of galaxies 

1 

0.1<z<0.5 

9.14 

9338 

2 

0.5<z<0.8 

9.47 

11760 

3 

0.8<z<l.l 

9.70 

13885 

4 

l.l<z<1.5 

9.93 

13640 

5 

1.5<z<2.0 

9.97 

12217 

6 

2.0<z<3.1 

9.97 

12641 
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Fig. 7.— NUV-r + versus r + -J color-color plots used to select 
quiescent and star-forming populations in our mass complete sam¬ 
ples at z=0.1-3.1. Galaxies with their rest-frame color NUV-r + > 
3.1 and NUV-r+ > 3(r+-J)+l are selected as quiescent systems 
Hilbert et al.ll20l 3f). 

The mass completeness limits for all the galaxies (quies¬ 
cent & star-forming) and quiescent systems are shown in 
Figure [6] As mentioned already, only less than 5% of 
galaxies could be missed in the lower mass population of 
galaxies. The mass completeness limit is higher for quies¬ 
cent systems. Here, we define six mass complete samples 
out to z^3, as shown in Figure[6] The properties of these 
mass complete samples are given in table 5. 

7.1.2. Selection of Quiescent/Star-forming Systems 

The selection of quiescent and star-forming systems 
can be performed using a single rest-frame color through 
the Color-Magnitude Diagram (CMD). However, the sin¬ 
gle color selection is problematic for several reasons. The 
existence of dusty, star-forming galaxies that mimic the 
color of quiescent systems can significantly contaminate 
the quiescent population and tend to unrealistically in¬ 
crease the qu iescent fraction and t heir comoving num¬ 
ber densities dBrammer et al.ll2009l) . The larger scatter 
in the rest-frame color at higher redshifts (mostly due 
to photo-z uncertainties) can wash out the red sequence 
and the disap pearance o f the rest- frame color bimodal- 
ity at z>1.5 (Will iams et al.l [200911 . Here, we use rest- 
frame two-color NUV-r + versus r + -J in order to select 
quiescent and star-forming populations in our mass com¬ 
plete samples. It has been shown that the rest-frame 


NUV-r + color is a better i ndicator o f the recent star 
formation activity dMartin et al.l (2007 ) and has a wider 
dynamical range com pared to t he more commonly used 
rest-frame U-V color (lllbert et al.ll201 3). Here, galaxies 
with their rest-frame color NUV-r + >3.1 an d NUV-r + > 
3(r + - J)+l are selected as quiescent systems (lllbert et al.l 
I2013T) . Figure [7] shows the rest-frame NUV-r + versus r + - 
J distribution of galaxies in the mass complete samples 
and the color cuts used to separate quiescent and star¬ 
forming systems. Two distinct populations of galaxies 
are clearly seen out to z~3. We stress that in our two- 
color selection technique, adding dust to the star-forming 
galaxies causes them to move diagonally from the bot¬ 
tom left to the top right of Figure [3 makin g them sep¬ 
arable from the q uiescent systems ((Williams et al.l[20091 : 
lllbert et al.ll2013j l. 


7.2. Evolution of Rest-frame Color as a Function of 
Environment For Quiescent Galaxies 

We now investigate the effect of environment on the 
rest-frame u*-Bj color of quiescent systems with differ¬ 
ent stellar masses and its evolution with redshift. We 
note that in selecting the quiescent population, we do 
not use the single rest-fra me u* -B./ color due to the is¬ 
sues expressed in Section 17.1.21 The selection is alter¬ 
natively done based on rest-frame NUV-r + versus r + -J 
plot. Here, we define two environments: high-density en¬ 
vironment (galaxy group and cluster scales) is defined as 
regions with overdensity values log(l+i5)>0.5 and low- 
density environment (field-like environment) is defined 
as those with log(l+<5)<0.5. The selection of the cut 
between low and high-density environments is somewhat 
arbitrary. However, it is shown (Darvish et al. 2015 in 
prep.) that the environmental effects (e.g., the increase 
in the fraction of quiescent galaxies with overdensity) 
start to effectively rise at log(l+<5)>0.5 in the COSMOS 
field. Throughout this work, we use this cut to separate 
the density field into low- and high-density environments. 
Based on this definition, we study the redshift evolution 
of the rest-frame u*-B,/ color as a function of environ¬ 
ment and stellar mass. Figure [8] shows the redshift evo¬ 
lution of the median rest-frame u*-B j color for quiescent 
galaxies with different stellar masses located in high- and 
low-density environments. For clarity, the evolution is 
shown for only a few stellar mass bins (AM=±0.1 dex 
around the selected masses). We choose the center of the 
redshift bins (table 5) as the redshift of the given points. 
The median color at each given environment, stellar mass 
and redshift is estimated using all the quiescent systems 
that are in that environment, with their redshift located 
in the redshift range of the mass complete samples and 
with stellar masses within AM±0.1 dex of the given stel¬ 
lar mass. The color uncertainties are estimated using 
10000 bootstrap resamples added in quadrature to the 
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Fig. 8. — Redshift evolution of the median rest-frame u*-Bj color for quiescent galaxies with different stellar masses located in high 
and low-density environments. For clarity, the evolution is shown for only a few stellar masses. The color uncertainties are estimated 
using 10000 bootstrap resamples added in quadrature to the median observational uncertainties in color. We assume that the evolution 
of rest-frame u*-Bj color for quiescent systems is linear with redshift and fit a straight line to the median colors at any given mass and 
environment. At a fixed redshift and stellar mass, the color of quiescent galaxies is independent of environment. However, at a fixed 
redshift and environment, the color of quiescent galaxies depends on stellar mass. Quiescent galaxies become redder with cosmic time and 
their evolution is independent of the environment they reside in. Quiescent galaxies more massive than log(M/M 0 )>lO become ~0.2 mag 
redder in rest-frame u*-Bj color since z~2.5. Since z~l, less massive systems (log(M/MQ)~9.5-10) redden by ~0.1 mag, regardless of 
their environment. We also find that more massive quiescent galaxies (log(M/M 0 )~ll) are ~0.1 mag redder in the rest-frame u*-Bj color 
compared to less massive (log(M/M 0 )~lO) systems and this color difference is almost independent of the environment and redshift. 


median observational uncertainties in color. We assume 
that the evolution of the rest-frame u*-B j color for quies¬ 
cent systems is linear with redshift and fit a straight line 
to the median colors at any given mass and environment: 

(u* - Bj)rest{z) = az + (u* - Bj) rest (z = 0) (12) 

The result of the fit to the median colors for different en¬ 
vironments and stellar masses is given in table 6. Accord¬ 
ing to Figure [5] and table 6, quiescent galaxies become 
redder with cosmic time and their evolution is indepen¬ 
dent of the environment they reside in. We particularly 
find that irrespective of environment, quiescent galax¬ 
ies more massive than log(M/M 0 )>lO become ^0.2 mag 
redder since z~2.5. This is in agreement with lKriek et all 
( 2008 ) who showed a reddening of ^0.25 mag in the rest- 
frame U-B color for massive (2xl0 n Mq) red galaxies 
in the field since z~2.3. This reddening is also seen in 
less massive quiescent systems (log(M/M e )=10-ll) and 
is independent of environment since z^3. Due to the in¬ 
completeness in stellar mass, we can not study the color 
evolution of less massive galaxies (log(M/M 0 )<lO) at 


z>l. However, at z<l, these systems (log(M/M 0 )~9.5- 
10 ) redden by ~0.1 mag, regardless of their environment. 
We also see that at any given redshift out to z~3, more 
massive quiescent systems are redder compared to less 
massive systems, regardless of their environment. Par¬ 
ticularly, we find that more massive quiescent galaxies 
(log(M/M 0 )~ll) are ~0.1 mag redder in the rest-frame 
u*-Bj color compared to less massive (log(M/M 0 )~lO) 
systems and this color difference is independent of the 
environment. The fact that more massive quiescent sys¬ 
tems are redder than the less-massive ones is a mani¬ 
festation of the Color-Magnitude Relation (CMR) since 
the stellar mass is proportional to the luminosity (mag¬ 
nitude) of galaxies. The CMR is seen in both the field 
(less dense) and cluster (dense) quiescent galaxies, and in 
the l ocal universe out to higher redshifts (see for exam- 


pie, ( 

3ower et all 1992: Blakeslee et al. 2003t Bell et al. 

2004: 

Mei et al. 2009; 

Wilson et al. 2009: Brammer et al. 

2009; 

Papovich et al.l 

2010 )). 


For quiescent galaxies, we find that the rest-frame color 
at a fixed environment depends on stellar mass but at a 
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Table 6: Parameters of the linear fit 
((u*-Bj)rest(z)=o:z+(u*-Bj) r est (z=0)) of the rest-frame u*-Bj color 
evolution for quiescent galaxies with different stellar masses. 


M 

log(M 0 ) 

OL 


(u*-B j) 

rest (z — 0) 

high-density 

low-density 

high-density 

low-density 

11.3 

-0.089±0.020 

-0.076±0.010 

0.969=1=0.016 

0.9564=0.012 

10.8 

-0.082±0.032 

-0.076±0.011 

0.917=1=0.026 

0.9204=0.016 

10.5 

-0.082±0.023 

-0.065±0.016 

0.916=1=0.024 

0.8824=0.019 

10.1 

-0.089±0.020 

-0.095±0.022 

0.882=1=0.018 

0.9014=0.018 

9.8 

-0.178±0.041 

-0.151=1=0.035 

0.8964=0.025 

0.8634=0.024 

9.6 

-0.176 

-0.157 

0.885 

0.853 


fixed mass is independent of environment. At first glance, 
this seems to contradict the studies showing a tight cor¬ 
relation between color and l ocal density of galaxies (for 
example. iBalogh et al.112004 ). However, the stellar mass 
also depends strongly on density and color-density rela¬ 
tion is actually a manifestation of a more fundamental 
color-mass relation. When controlled for the stellar mass 
(or luminosity), the color-density relation becomes inde¬ 
pendent of the environment (overdensity). 

The environmental independence of color-mass relation 
is consistent with several p revious s tudies. In the local 
universe, For example, IBalogh et al.1 ( 2004 .) showed that 
at a fixed luminosity (equivalent to a fixed stellar mass), 
the mean rest-frame color of red galaxies in the SDSS at 
z<0.08 is nearly i ndepend e nt o f their environment. In¬ 
dependently, [Baldrx^LaLj ( 200 6) showed that at z<0.1, 
the color-mass distribution of galaxies (red & blue) does 
not significantly change with environment. The same 
result is also seen in iHogg et all (l2004tl who found the 
independence of red bulge-dominated galaxy colors on 
environ ment in the SDSS at z<0.12. Using the SDSS 
data, iPark et ahl ( 2007 ) showed that when morphology 
and luminosity (equivalent to stellar mass in our study) 
are fixed, other physical properties of galaxies, including 
color is nearly independent of local density of galaxies. 
They also showed that more luminous galaxies are red¬ 
der than less luminous systems and this is independent 
of the environment. This is consistent with our results. 
The environmental invaria nce o f color for bright galaxies 
in the SDSS is also seen in lTanaka et al.1 (120041 ). 

Several studies at higher redshif ts are also consistent 
with our results. For example, iCassata et al.1 (20071 ) 
showed that at z^0.7 in the COSMOS field, the ob¬ 
served color-magnitude diagram for red and early-type 
galaxies is independen t of the local density of galaxies. 
iScodeggio et al. (200ft ) found that for a sample of galax¬ 
ies in the VVDS survey at 0.2<z<1.4, color is indepen¬ 
dent of the local density of galaxies at a fix ed st ell ar mass. 
Using the zCOSMOS data out to z~l, iMoresco et al.1 
(2010) found that the rest-frame color of red galaxies at 
a fixed stellar mass is almost independent of environ¬ 
ment but at a fixed environment, depends on the stellar 
mass. They also found that the average colors of mas¬ 
sive red galaxies (log(M/M 0 )=lO.8) are redder than low- 
mass galaxies ((log(M/M 0 )=lO) throughout their entire 
redshift range. This is entirely consistent with our re¬ 
sults in this section. By extending this study to higher 
re dshifts. we d emonstrate that the results hold to z~3. 
iCucciati et al.l (l2010h showed that for red massive galax¬ 
ies ((log(M/M 0 )>lO.7) in the zCOSMOS lOk-sample, 
the color-density relation is globally flat up to z~l, con - 
sistent with our results. Recently, iBassett et al.l (20131 ) 


compared the rest-frame U-B color of galaxies in a cluster 
with those in the field at z=1.6 in the CANDELS-UDS 
field and found no difference between the color of quies¬ 
cent galaxies in these two environments. This agrees well 
with our results at z>l. 

7.3. Effect of the Environment on the Comoving 
Number and Mass density of Massive Galaxies 

In this section, we investigate the evolution of comov¬ 
ing number (n) and mass ( p ) density for massive (>10 n 
M 0 ) galaxies (quiescent & star-forming) in different en¬ 
vironments since z~3. We start our analysis by study¬ 
ing the evolution of n & p for massive (> 10 11 M 0 ) 
quiescent/star-forming system regardless of their envi¬ 
ronment. This allows us to compare n & p with previous 
studies which often do not consider the role of environ¬ 
ment in their analysis. Later, we discuss the role of en¬ 
vironment. 

Figure [9] shows the evolution of comoving number and 
mass density of massive (> lO n M 0 ) quiescent and star¬ 
forming galaxies since z~3, along with comparison with 
some previous studies. Some of these studies use Kroupa 
IMF to estimate stellar masses. For those, we modify 
the mass densities based on the Chabrier IMF (M chab 
~ 0.89 M Kroupa)- The redshift of our data points is se¬ 
lected as the center of redshift bins introduced for mass 
complete samples (Table 5). Error bars for number den¬ 
sities incorporate both Poisson error and uncertainties 
due to the cosmic variance. Uncertainty in mass density 
is estimated using 10000 bootstrap resamples, added in 
quadrature to the observational uncertainties and those 
due to the cosmic variance. For observational uncertain¬ 
ties, we use the median of the stellar mass uncertainties of 
the galaxies in each su bsamp le. We use the cosmic vari¬ 
ance calculator of iMoster et all (20111 ) to estimate the 
cosmic variance uncertainties. For ll<log(M/M 0 )<11.5 
galaxies, the fractional uncertainties due to the cosmic 
variance change between ~15-10% at ~z=0.1-3 in the 
COSMOS. 


7.3.1. Massive Quiescent Galaxies 

Figure 0 (a) and (b) show a rapid increase in the co¬ 
moving number and mass density of massive quiescent 
systems from z~3 to z~l. However, we observe a change 
of pattern at z~l. From z~l to the present time, the 
number & mass density of massive quiescent galaxies re¬ 
main almost unchanged within uncertainties. From z=2- 
3.1 to z=l.5-2.0, the number and mass density of massive 
quiescent galaxies are increased by a factor of ~2.5. From 
z=l.5-2.0 to z=0.8-1.1, the rise in the number and mass 
density of massive quiescent galaxies is steeper and they 
are increased by a factor of ~6-8. 

The lack of evolution (or insignificant evolution) since 
z~l in the comoving number and mass density of massive 
quiescent systems is consistent with several studies. In 
terms of sa mple selection the clo sest works to our stud y 
are those of lllber t et al.1 (201 3D fc iMuzzin et al.l (2013bD . 
Illbert, etHl (|2013l) did not find any significant evolution 
in the number and mass density of high-mass end qui¬ 
escent galaxies at z<l, cons istent wi th o u r resu lts. As 
seen in Figure [9] (a) & (b), IMuzzin et al.1 (l2013bl) data 
points follow our observed trends for n & p relatively 
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Fig. 9.— Evolution of comoving number (n) and mass (p) density for massive (>10 11 Mg) galaxies (quiescent & star-forming) since z~3, 
along with comparison with some previous studies. Error bars for number densities incorporate both Poisson error and uncertainties due to 
the cosmic variance. Uncertainty in mass density is estimated using 10000 bootstrap resamples, added in quadrature to the observational 
uncertainties and those due to the cosmic variance. The comoving number and mass density of massive quiescent systems rapidly increase 
from z~3 to z~l. However, from z~l to the present time, they remain almost unchanged within uncertainties, we find almost no evolution 
in the comoving number and mass density of massive star-forming galaxies since z~2. The number & mass density of massive star-forming 
galaxies are slightly increased from z=2.0-3.1 to z=l.5-2.0 (a factor of ~2) and remain almost unaltered since then. 


TABLE 2 

Table 7: Comoving number density, n (in units of 10 -4 Mpc -3 ) & mass density, p (in units of 10 7 M@ Mpc -3 ) of massive 

(> 10 11 Mg) QUIESCENT & STAR-FORMING GALAXIES IN DIFFERENT REDSHIFT BINS. 


Property 



redshift 



0.1<z<0.5 

0.5<z<0.8 

0.8<z<l.l 

l.l<z<1.5 

1.5<z<2.0 

2.0<z<3.1 

n Q 

2.17T0.37 

2.08T0.30 

2.68T0.35 

1.05T0.13 

0.37T0.05 

0.15T0.02 

n SF 

1.22±0.25 

0.97T0.15 

0.84±0.12 

0.91T0.11 

0.88T0.11 

0.44±0.05 

PQ 

4.20T1.17 

3.73T1.00 

4.57T1.20 

1.71T0.44 

0.61T0.16 

0.23T0.06 

PSF 

1.78T0.55 

1.40T0.37 

1.32±0.35 

1.37T0.35 

1.31T0.34 

0.72±0.18 


well, within uncertainties. As seen in Figure [9] (a) & 
(b), for quiescent galaxies at z<l, our result is also in a 
good agreement with >1O 1:L M 0 red-se qu ence galaxies in 
C0MB017+4 Survey from iNicol et all (|2011l) , massive 
U-V-J selected quiesc ent galaxies in NMBS survey from 
iBrammer et al.l (|20:H), log(M/MM ) =ll-1 1.5 red galaxies 
in zCOSMOS from iPozzetti et a ll (I2010D , quiescent sys - 
tems in PRIMUS survey from~ lMoustakas et al.l ( 1201 30 
and >1O U M 0 rest-frame (U-V ) selec ted red galaxies in 
VIPERS from iDavidzon et, al.l (120131) . Several luminos ¬ 
it y function studies a t z<l such as lScarlata et, al.1 (|2007t ) 
& iCool et all (120081) have seen almost no evolution in 


the number density of very luminous (L~2.5-3 L*) red 
or early-type galaxies since z~l, in agreement with our 
results. 

The sharp rise in the comoving number and mass den¬ 
sity of massive quiescent systems from z=3-l is also 
in a greement with som e previous studies. For exam¬ 
ple, illbert et~ahl (j2013f ) found that the comoving num¬ 
ber and mass density of 10 11 M 0 quiescent galax¬ 
ies increase by factors 25 & 13, respectively, between 
2.5<z<3 & 0.8<z<l.l. This is similar to our result 
showing an increase of ^15-20 in n & p from z~3 
to z^l for massive quiescent systems. Our results 
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2011 


2011 


also agree qualitatively/quan t itatively with dNicol et al 
IBrammer et all 1201 lk iDominguez Sanchez et al 
iMuzzin et alJl2013bfl 


7.3.2. Massive Star-forming Galaxies 

The situation is different for massive (>10 n M 0 ) star- 
forming systems. According to Figure [9] (c) and (d), 
we find almost no evolution in the comoving number 
and mass density of star-forming galaxies since z~2. 
The number & mass density of massive star-forming 
galaxies are slightly increased from z=2.0-3.1 to z=1.5- 
2.0 (a factor of ~2) and remain almost unaltered since 
then . This sli ght (i nsignificant) increase is seen in 
IBrammer et al.l (|2011f ) from z~2.5-2.0, followed by de¬ 
cline since z~2. We argue that we do not find any sign 
of decline in the comoving number & mass density of 
massive star-forming systems since z^2 in our st udvU 


Muzzi r et al. 2013U Ilbert et al. 2013 

Davidzon et al. 

2013: Sobral et al.l I2014H. However, 

Moustakas et al. 


(1201311 found that the number (& mass) density of mas¬ 


sive star-forming galaxies are decreased slightly by a fac¬ 
tor ~2 since z~l (Figure|9](c) & (d)). Within uncertain¬ 
ties, we do not find a significant disagreemen t between 
our result and that of lMoustakas et all 1 2013 ). although 
we detect no evolution in n & p for massive star-forming 
galaxies. We argue that even in the presence of a real 
evolution in number and mass density of massive star¬ 
forming systems since z~l, we would not be able to ob- 
se rve it due to the la r ger er ror bars in our study compared 
to lMoustakas et ahl ( 2013( 1 . 


The comoving number and mass density values for mas¬ 
sive (> 10 11 M 0 ) quiescent/star-forming populations are 
given in table 7. 


7.3.3. Massive Quiescent/Star-forming Systems in Different 
Environments 

Now, we investigate the role of environment in shaping 
the number and mass density of massive quiescent and 
star-forming galaxies. In order to estimate comoving 
number and mass densities in different environments, 
one needs to determine what fraction of the density 
field is occupied by massive galaxies in different envi¬ 
ronments. This is not practically straightforward as the 
definition of environment and the selection of low- and 
high-density regions are arbitrary. We highlight that we 
only wish to relatively compare the number and mass 
densities in different environments and not to determine 
the exact values in these regions. If there are any 
systematic errors in determining the number and mass 
densities in different environments, they would likely be 
cancelled out as we take the ratio of these in low- and 
high-density regions. We define the ratio of comoving 
number density of massive (> 10 11 M 0 ) galaxies in 
high-density (n^) to low-density (n id) environments as: 

n hd Ei i/v; M 

nw Zjl/Vj* [ ’ 

3 IBrammer et al.l (120111) state in their paper that the comoving 
number density evolution of >10 11 Mq star-forming galaxies is 
nearly flat out to z=2.0 


where and V ld are the volumes associated with the 
ith and jth massive galaxy in high- and low-density 
environments, respectively. The volume assigned to 
each galaxy is estimated using its Voronoi area (A) and 
the radial comoving length (Al) that corresponds to its 
z-slice; i.e. Vj=AiAb. Similarly, we define the ratio 
of mass density of massive (> 10 11 M 0 ) galaxies in 
high-density (phd) to low-density (pu) environments as: 

Phd = E,M dd /V dd 

pid EjM l //v‘ d [ 

where M; and M ld are the stellar masses associated with 

1 J 

the Zth and jth massive galaxy in high- and low-density 
environments, respectively. 

Figure [TO] shows the redshift evolution of n^d/n u and 
Phd/Pid for massive quiescent and star-forming galaxies. 
Uncertainties in number density ratios incorporate both 
Poisson error and the cosmic variance. The error bars for 
mass density ratios are estimated using 10000 bootstrap 
resamples, added in quadrature to the observational un¬ 
certainties in stellar mass and those due to the cosmic 
variance. For massive (> 10 11 M 0 ) star-forming sys¬ 
tems, the number and mass densities in different envi¬ 
ronments remain almost the same and we find almost 
no evolution in comoving number and mass density ra¬ 
tios in different environments with cosmic time. Massive 
star-forming galaxies populate dense and less-dense re¬ 
gions almost equally, regardless of their redshift. The 
situation is different for massive quiescent galaxies. At 
z>1.3, there is no significant evolution in the comov¬ 
ing number and mass density ratios in different environ¬ 
ments. At these redshifts, the number and mass density 
ratios for massive quiescent galaxies are almost equal to 
those of massive star-forming systems. However, at lower 
redshifts (z<1.3), we find a significant evolution in the 
comoving number and mass density ratios for massive 
quiescent systems in different environments. These ra¬ 
tios for massive quiescent galaxies monotonically increase 
with cosmic time at z<1.3. At z<0.5, the number, as 
well as mass density of massive quiescent galaxies are 
~1 dex higher in denser regions compared to less dense 
environments. Dense environments at low redshifts are 
populated by massive quiescent galaxies. 

We stress that part of the evolution in the comoving 
number and mass density of massive quiescent systems 
in different environments might be due to the growth of 
the large scale structure with cosmic time. As time pro¬ 
gresses, less dense environments eventually coalesce to 
assemble more massive, denser regions. We also mention 
that the lack of evolution at higher redshifts in the num¬ 
ber and mass density ratios for massive quiescent galaxies 
in different environments does not necessarily mean that 
the environment at higher redshifts is not able to quench 
the star formation activity in galaxies as effectively as it 
does at lower redshifts. The likelihood of finding massive 
halos (that are able to suppress the star formation activ¬ 
ity) at higher redshifts is low, given the small size of the 
COSMOS field and the rarity of massive halos at higher 
redshifts. If we could numerously find such dense mas¬ 
sive halos at higher redshifts in our survey, we might be 
able to see similar environmental trends that we observe 
in lower redshifts and the local universe. 
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(a) 


(b) 




Fig. 10.— (a) The ratio of the comoving number density in high-density environments to low-density environments (n^/n id) for massive 
(> 10 11 Mg) quiescent and star-forming galaxies as a function of redshift. (b) The ratio of the stellar mass density in high-density 
environments to low-density environments ( phd/pld .) for massive (> 10 11 Mg) quiescent and star-forming galaxies as a function of redshift. 
Both number and mass density ratios do not significantly change with redshift for massive star-forming galaxies. For massive quiescent 
systems, these ratios significantly change at lower redshifts, indicating the prevalence of massive quiescent galaxies in denser environments 
at lower redshifts. 


It is challenging to make a direct, quantitative compari¬ 
son between our results in this section and similar stud¬ 
ies in the literature, due to the differences in sample se¬ 
lections, the definition of environment and the classifi¬ 
cation of low- and high-density regions. However, we 
qualitatively compare our results in this section with 
some similar studies. Using zCOSMOS data at z<l, 
iBolzonella et al.l 1 2010 ) studied the galaxy stellar mass 
function (SMF) as a function of galaxy type and envi¬ 
ronment. They found that massive galaxies are pref¬ 
erentially located in high-density environments, charac¬ 
terised on average by a higher M* and that the spectral- 
type selected early-type galaxies dominate the high-mass 
(log(M/MQ)>10.7) end of the SMF. In other words, mas¬ 
sive,red galaxies are preferentially reside in high-density 
regions at z<l, in agreement with our results for mas¬ 
sive, qui escent g alaxies in dense environments. Recently, 
iMortlock et al.l (|2015l) studied the SMF of galaxies in 
UKIDSS/UDS and CANDELS fields and found that at 
higher redshifts (z>l), the SMF parameters for galaxies 
in low and high densities are almost the same. How¬ 
ever, at z<l, they found evidence that: 1) the high-mass 
end of the galaxy SMF is more dominated by galaxies in 
dense environments. 2) The high-mass end of the Sersic 
index > 2.5 SMF is dominated by quiescent galaxies. 3) 
on average, </>* and M* values of the SMF are larger for 
red galaxies in denser regions compared to those in less- 
dense environments (At z~0.5, cf>* is ~1 dex larger for 
red galaxies in denser regions, similar to our result in 
this section). Combining all these together, we conclude 
that at lower redshifts, high-density environments cause 
the build-up of high-mass, quiescent galaxies. This is 
completely consistent with our results in this section re¬ 
garding the higher number and mass densities in denser 
regions for massive, quiescent systems at lower redshifts 
and the lack of any significant differences at higher red¬ 
shifts. 

We stress that our results here rely on the density field 
estimation based on the flux-limited sample (K s <24) in¬ 
troduced i n Sect ion [3 However, as we already discussed 
in Section 17.1.11 this magnitude cut results in a sam¬ 


ple with varying stellar mass as a function of redshift. 
In other words, the “typical” galaxies defining the envi¬ 
ronment change with redshift, with less massive galaxies 
defining the low-z and more massive systems defining 
the high-z environments. Since there is some degree of 
correlation between galaxy stellar mass and environment 
I Kauffmann et al.il 2004 iBaldrv et all I2006D , our results 
in this section might be affected by selection of galaxies 
used to estimate the density field (environment). How¬ 
ever, in the Appendix, we investigate the role of different 
samples used for density estimation, on the results and 
will show that the results do not change (even become 
stronger) and we recover the same trends discussed in 
this section. 

We conclude that the comoving number and mass density 
of massive, star-forming systems do not evolve much with 
redshift, regardless of their environment. This scenario 
is also true for massive, quiescent galaxies at higher red¬ 
shifts (z>1.3). However, at lower redshifts, the comoving 
number and mass density of massive, quiescent galax¬ 
ies are greater in high-density environments compared 
to less-dense regions. This highlights the significant role 
of environment in quenching the star formation activity 
in galaxies at lower redshifts. 

8. SUMMARY & CONCLUSIONS 

In this work, we used a K s -band selected sample of 
galaxies with accurate photometric redshifts in COS¬ 
MOS at z=0.1-3.1 in order to estimate the density field. 
The density field was determined with the weighted 
versions of the adaptive kernel smoothing, 10 t?l & 5 th 
NN, Voronoi tessellation and Delaunay triangulation 
methods. We evaluated the performance of each density 
estimator using extensive realistic and Monte-Carlo 
simulations. We later defined two environments and 
studied the effects of environment on mass complete 
sample of quiescent and star-forming galaxies out to 
z~3. The Rest-frame NUV-r + -J color-color plots were 
used to separate the galaxy population into quiescent 
and star-forming systems. We investigated the redshift 
evolution of the rest-frame u*- B j color for quiescent 
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galaxies as a function of environment and stellar mass. 
We also studied the evolution of comoving number and 
mass density of massive quiescent and star-forming 
galaxies and their dependence on the environment of 
galaxies. Our main results are: 

1- We find an overall good agreement between the 
density field estimated with weighted versions of adap¬ 
tive kernel smoothing, nearest neighbor (NN), Voronoi 
tessellation and Delaunay triangulation methods over 
~2 orders of magnitude. 

2- Extensive simulations show that the adaptive kernel 
smoothing and Voronoi tessellation outperform other 
methods in estimating the density field of galaxies. We 
recommend using these estimators as a more reliable 
and robust substitute for the widely-used NN or count 
in aperture methods. 

3- At fixed stellar mass, the median rest-frame u*-Bj 
color of quiescent galaxies is independent of the environ¬ 
ment they reside in. Quiescent galaxies become redder 
with cosmic time and their color evolution is indepen¬ 
dent of their environment. Since z^3, more massive 
quiescent galaxies (log(M/Mg)>10) have become ~0.2 
mag redder in rest-frame u*-Bj whereas less massive 
quiescent systems have reddened by ~0.1 mag since z~l. 

4- On average, more massive quiescent galaxies are 
redder compared to less massive ones at any given 
redshift, regardless of their environment. The lack 
of a correlation between color and environment at 
fixed stellar mass for quiescent galaxies suggests that 
the relation between stellar mass and local density of 
galaxies is more fundamental than the color-density 
relation. 

5- The average comoving number and mass density of 
massive (log(M/M Q )>ll) star-forming galaxies have not 
evolved much since z^3. However, for massive quiescent 
galaxies, number and mass densities sharply rise from 
z~3 to z~l and remain almost unchanged since then. 

6 - The evolution of comoving number and mass density 
of massive star-forming galaxies do not depend on their 
environment. They remain almost unchanged since z~3, 
regardless of their host environment. The situation is 
different for massive quiescent galaxies. The comoving 
number and mass density of massive quiescent galaxies 
do not change much with environment and redshift from 
z~3 to z~1.3, similar to those of star-forming galaxies. 
However, at lower redshifts (z<1.3), we find a significant 
evolution in the number and mass density of massive 
quiescent galaxies in denser environments compared 


to less-dense regions. Dense environments at lower 
redshifts are populated with massive quiescent galaxies. 
This signifies the role of environment in quenching the 
star formation activity in galaxies at lower redshifts. 

This paper is the first one in a series and provided the 
required tools for the density-based environmental study 
of galaxies. In a following paper (in prep.), we will study 
the effects of local environment of galaxies on their SFR, 
sSFR, rest-frame color and the fraction of quiescent/star¬ 
forming systems as a function of redshift, stellar mass 
and galaxy type. We will also discuss the fractional role 
of stellar mass and environment in suppressing the star 
formation activity in galaxies. 
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APPENDIX 

EFFECT OF DENSITY FIELD ESTIMATION BASED ON DIFFERENT SAMPLE SELECTIONS 

In Section [2j we used a flux-limited (K s <24) sample of galaxies in order to estimate the density field. All the 
main results presented in Sections 17.21 & 17.3.31 were based on this flux-limited sample. Here, we try two other samples 
of galaxies to estimate the density field and investigate their effects on the results presented in Section 17.3.31 Same 
as before, we use the Voronoi tessellation method as the density estimator and the definition of low- & high-density 
environments given in Section [7.21 These new samples are defined as follows: 

1. In addition to conditions given in Section [2l we select all galaxies that are more massive than the stellar mass 
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(a) (b) 




Fig. 11. — Similar to Figure llOl but the results are shown for different samples used to estimate the density field of g alax ies. Filled 
circles, empty squares and empty triangles show the trends based on a flux-limited sample (K s <24, same results in Figure ITOt) . a sample 
of galaxies with (K s <24 Sz M>M comp ) and a sample of galaxies with ( K s <2 4 Sz log(M/M0)>9.97, similar to a volume-limited sample), 
respectively (see Appendix). We see the same trends obtained in Section f7.3.3l Our results are not significantly influenced by the selection 
of the samples used to estimate the density field of galaxies. 

completeness limit of samples defined in Section f7. 1.II and table 5 (K, <24 & M>M comp ). 

2. In addition to conditions given in Section O we select all galaxies that are more massive than the stellar mass 
completeness limit of highest redshift sample defined in Section !?. 1. H and table 5 (K s <24 & log(M/M 0 )>9.97). 
This is similar to a volume-limited sample of galaxies. 


Using these new samples, we re-estimate the density fields and re-investigate the results in Section 17.3.31 Figure ITT] 
is similar to Figure [TD1 but it also shows the results based on the new samples us ed to estimate the density field. 
According to Figure fill we are able to retrieve the main trends obtained in Section [7.3.31 for quiescent (Q) and star¬ 
forming (SF) galaxies. However, we mention that for massive quiescent galaxies, the number and mass density trends 
become even more amplified at z<l when we use the new samples for density estimation. This is expected since 
the new samples target more massive galaxies in order to estimate the density field, which tend to be more strongly 
clustered compared to less massive systems. We highlight that our results in Section [7.3.31 are not significantly affected 
by the selection of the samples we used for density field estimation. 
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